function inv_hs  = inv_hs(shares, rs)
    global elasts
    inv_hs = ((shares ).^(1./elasts) )./rs ; 
    
end